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Abstract 

A bottleneck simulation of road traffic on a loop, using the deterministic cel- 
lular automata (CA) Nagel-Schreckenberg model with zero dawdling prob- 
ability, reveals three types of stationary wave solutions. They consist of i) 
two shock waves, one at each bottleneck boundary, ii) one shock wave at the 
boundary and one on the "open" road, and hi) the trivial solution, i.e. ho- 
mogeneous, uniform flow. These solutions are selected dynamically from a 
range of kinematicly permissible stationary shocks. This is similar in fash- 
ion to the wave selection in a bottleneck simulation of the optimal-velocity 
(OV) model, which is explained by a travelling wave phase-plane analysis 
of the corresponding continuum model. It is yet another strong indication 
that CA and OV models share certain underlying dynamics, although the 
former are discrete in space and time while the latter are continuous. 

Keywords: Traffic flow, bottleneck, cellular automata, Nagel-Schreckenberg 
model, optimal velocity, pattern formation 
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1 Introduction 



Cellular automata (CA) models have been widely used to simulate traffic 
flow on highways and road networks [HElEllllElElEllHlE], in particular the 
Nagel-Schreckenberg model [101 UT] . Together with car- following (CF) and 
continuum models, they represent the three predominant classes of traffic 
models. 

Analytical work by Berg et al [12] and Lee et al [13] has established a 
link between car-following models based on ordinary differential equations, 
and continuum models based on partial differential equations. While an 
analytical link between CA models and either CF or continuum models 
is still missing (mean- field theory aside), the dynamics of all three classes 
exhibit many common features such as sub-critical bifurcations, limit cycles 
and pattern formation [T4l [T5] . 

Bottlenecks are the major cause for highway congestion and, therefore, 
have been studied in some detail [TH [T71 HH EE EDI EI]- I n this paper, a 
wave selection analysis of a bottleneck simulation reveals a fundamental link 
between the dynamics of CA models and optimal- velocity (OV) models [22] . 
which belong to the class of CF models. 

The paper is organised as follows. We begin with the description of the 
CA bottleneck simulation in section El Section [3] presents the resulting wave 
patterns and how they can be interpreted in the fundamental diagram. In 
the thermodynamic limit, boundaries between the different wave regimes are 
computed. Section [4] briefly addresses the impact of the dawdling probability 
on the effective fundamental diagram. 

Section [5] discusses a bottleneck simulation in the car-following OV model, 
resulting in identical wave patterns as compared to the CA model. In con- 
trast to the CA model, analytical tools can be applied to explain the wave 
selections. In particular, continuum theory is used to investigate travelling 
waves in terms of plateau connections between equilibrium flow solutions. 
A strong link between CA and OV models emerges. 

Finally, some of the results are compared to work in the literature which 
is related to bottleneck simulations (section Ej) , before we draw some con- 
clusions and mention planned future work in the last section. 

2 Cellular Automata bottleneck simulation 

For traffic on a loop (periodic boundary conditions) of length L, a bottleneck 
of length Lb is located at < x < Lb (Fig.[T]). The system is simulated 
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Figure 1: The bottleneck simulation assumes periodic boundary conditions 
(traffic on a loop). The "open" road (Lb < x < L) is modeled with the 
conventional Nagel-Schreckenberg model [10] . The bottleneck (0 < x < Lb) 
uses the same model except for a smaller maximum speed y™ ax . In the 
simulations, we set Lb = 200 and L = 1000. 



with the Nagel-Schreckenberg (NS) CA model |10j for vanishing randomness 
(p = 0) and a reduction in top speed from y max = 5 on the "open" road to 
y max — 3 m the bottleneck. All other model parameters remain the same. 
Initially, N cars are randomly distributed along the road and the system is 
updated up to t = 10 6 time steps. We set Lb = 200 and L = 1000 and study 
the emerging wave patterns. We also choose p = in order to avoid the jam 
formation in the NS model, which would interfere with the stationary wave 
patterns and complicate the analysis. 

Figure [2] is a qualitative plot of the equilibrium velocity function used 
in the update rule of the NS CA model for both the bottleneck and the 
open road. Here, h! = 1, y max = 5 and v^ ax — 3. Since the density p is 
continuous in this plot but the CA model is discrete in space, p must be 
interpreted in the thermodynamic limit as an average density taken over 
many cells. 

Note that the wave selection on a loop is fundamentally different from 
wave selection on an open road with different conditions at the upstream and 
downstream boundary, respectively [SJEJEj. On an open road, travelling 
waves can only interact once with stationary structures since they do not 
move around a loop. Therefore, an open road gives rise to new wave solutions 
as compared to the loop, determined by the boundary conditions El Ej. 
This is a very different setup from what is presented in this work. 
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Figure 2: The equilibrium velocity function of the CA model, velocity V as a 
function of the headway h (distance) between the cars. It is a piecewise linear 
function in the bottleneck and on the open road with an upper bound of y max 
on the open road and y™ ax in the bottleneck. The functions are partially 
identical which translates into partially identical fundamental diagrams (see 
Fig.[3j). Here, h! represents the minimum distance between the cars. 



3 Wave selection in the fundamental diagram 

We will use the fundamental diagram (FD), i.e. flux versus density, and 
kinematic wave theory to interpret our numerical results. Based on the 
equilibrium velocity function in Fig. [2] and using the equilibrium flow-density 
relation of uniform flow (q: flux, p: density) 

q = v(p)p, p= 1/h, (1) 

the FD is shown for both the bottleneck 

„ _ / 3 p ; < p < V4, , 9 , 



and the open road 



1/4 < p< 1 

< P < 1/6, , , 

1/6 < p < 1 1 J 

in Fig. 13 The two curves merge into the same function for p > 1/4. 

Generally speaking, we could expect as many as six stationary wave so- 
lutions for a bottleneck on a loop as t — > oo. In the FD, five of them are 
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Figure 3: Stationary shock waves visualized in the fundamental diagram: 
From six possible wave connections (the trivial uniform flow solution in 
Figs.[7|and[8]is not shown here), only three emerge as dynamical solutions, 
displayed in Figs. [5)13 



visualized as chords with zero gradient due to the requirement of vanish- 
ing wave speed, based on kinematic wave theory |20| . By kinematic waves 
we refer to density waves described by the first-order hyperbolic Lighthill- 
Whitham model of traffic flow [23J 



For stationary waves, we have p t = and the equation can be integrated to 
yield 



Therefore, stationary shock waves, whose profiles are given by p(x), exhibit 
constant flux along the road, represented by a horizontal straight line in 
the fundamenatal diagram. In Fig. [HI each chord is a straight line between 
points on the fundamental diagrams, one point representing the bottleneck 
and one representing the open-road FD. They connect plateaus between the 
bottleneck and the open road (case 2, 3 and 4) in case of one stationary shock 
wave at each bottleneck boundary. However, they can also entail a plateau 
connection on the open road as in cases 1 and 5. These five stationary wave 
patterns are shown in Fig. [4] in terms of density distribution along the loop. 
In addition, there is the trivial wave solution of homogeneous uniform flow. 
In principle, we could think of further wave patterns but we will restrict the 
analysis to the simplest cases featured here. 



Pt + q(p) x = pt + q\p)p x = 0. 



(4) 



q(p(x)) = const. 



(5) 
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(a) 

P=°f 



(b) 
: p<0.25 



(c) 



(d) 



(e) 




Figure 4: Stationary wave patterns of Fig. El While cases (a) and (b) occur 
in the simulations (cases (1) and (2) in Fig. [3]), as shown in Figs. [5] and El 
cases (c), (d), and (e) do not emerge. This is a very close analogy to the 
wave pattern selection of a bottleneck simulation with the optimal-velocity 
model (see section [5} [24| [25]. 



We found in our simulations that only three wave patterns are selected 
from this range of possible solutions. They consist of the following (p: av- 
erage density on the loop): 

• Case 1: 0.17 < p < 0.25 

Stationary wave pattern that connects two plateaus by one shock at 
the downstream boundary of the bottleneck and one classical (Lax) 
shock on the open road (see Figs.SK and[5j): The resulting bottleneck 
headway (distance between cars) is exactly at ds = 4 and, therefore, 
the bottleneck is at maximum flow. On the open road we find the 
headway to be near d n = 7, or exactly d Q = 20/3 on average. 

• Case 2: < p < 0.17 

Stationary wave pattern that connects two plateaus by a shock at the 
upstream and downstream boundary of the bottleneck, respectively 
(see Fig. [4}3 and [6]): In the bottleneck > 4 and on the open road 
d Q > 20/3. It shall be stressed that it takes a very long time for 
the system to reach steady-state due to the small interaction of cars 
on the open road. Hence, Fig. [6] should be considered as a transient, 
quasi-steady state. 

• Case 3: p > 0.25 

Trivial flow solution, i.e. homogeneous, uniform flow: Unless the aver- 
age headway d is close to an integer number, as is the case in Fig. El the 
individual headways d n oscillate around the average headway d, exhib- 
ited in Fig. [8l However, this is an effect solely due to the discretization 
of space, and the flow solution can still be considered uniform. In 
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particular, the average flow measured in Fig. El q = 0.55, equals ex- 
actly the uniform flow in the FD, as expected from an average density 
p — 0.45, namely q = 1 — p — 0.55. 



d=20/3 




1000 



Figure 5: Stationary wave pattern of bottleneck simulation with average 
density 0.17 < p < 0.25 (here: p = 0.20): Two shocks emerge, one at the 
downstream bottleneck boundary and one classical (Lax) shock on the open 
road. The flow in the bottleneck is at its maximum. 



This leads to three open questions: 

1. Why is there a transition between the structures at p 
p = 0.25? 



0.17 and 



2. What determines the location of the classical shock on the open road 
in case 1? 

3. Why do we not observe the other wave patterns? 

We will now elaborate on all three questions. 

The two headways in Fig. [6] are determined by the conservation of cars 
and by imposing zero wave speed (zero gradient of the chord in the FD), de- 
scribed by Eq. ([5]). This can be written as two equations with two unknowns, 
the bottleneck headway ds and the open-road headway d Q . Neglecting finite 
size effects, conservation of cars reads 



L — Lb ^_ Lb 
d d B 



N. 



(6) 
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Figure 6: Stationary wave pattern of bottleneck simulation with average 
density p < 0.17 (here: p = 0.142): Two shocks emerge, one at each bottle- 
neck boundary. 



We find for the wave speed criterion (equal fluxes qs and q Q in both road 
segments) 

3 5 

qB = q ^ 3pB = 5p =>> — = — , (7) 

dB d Q 

where pB and p Q denote the bottleneck and open-road density, respectively. 
The system ©-([7]) can be solved for ds- It yields 

d B = (8) 

In Fig.El we have L = 1000, L B = 200, N = 142 and, hence, d B = 4.79. 
This equals p B = 1/de = 0.21, which coincides with the numerical value. 
The value for d Q follows correspondingly. 

The maximum amount of vehicles that the wave structure in Fig. [6] can 
support, however, is reached when ds = 4.0 and, determined by zero wave 
speed, d Q = 20/3. For Lb — 200, this corresponds to an average density of 

This is coincidentally close to q = 1/6, the maximum of q Q , but varies with 
the choice of Lb- If we increase p beyond this, the wave pattern in Fig. [5] is 
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Figure 7: Trivial flow solution for bottleneck simulation with average density 
p > 0.25 (here: p = 0.333 d « 3). 



triggered, with a bottleneck headway exactly at ds = 4 and d G = 20/3 for 
the other plateau value. This is shown by chord 1 in Fig. 03 The length of 
the second plateau L p is now determined by the conservation of cars alone. 
Setting p = 0.2 in Fig.EJ we write 

P=\(l- L p /L) + ^(Lp/L) ^L p = 500, (10) 

which is very close to the numerical value of L w 510. Note that finite size 
effects impose limits on the accuracy of estimates. For p > 0.25, which is 
the maximum of the length L p of the open-road plateau equals zero and 
this wave pattern must vanish. We are left with the trivial flow solution of 
Figs. [3 and [8] since the average density exceeds the maximum density that 
can support the wave structure in Fig.0 

Thus, we have answered questions 1 and 2. Before we turn to question 3 
in section [5j let us study the effective fundamental diagrams obtained in the 
simulations. Here, we will only briefly discuss the impact of p > on the 
fundamental diagram and leave a more detailed study to future research. 

4 Impact of the dawdling probability on the 
fundamental diagram 

The well-defined wave patterns for p = overlap with instabilities of the flow 
when p > 0. Therefore, the dawdling probability p has a significant impact 
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Figure 8: Trivial flow solution for bottleneck simulation with average density 
p > 0.25 (here: p = 0.45 d « 2.22). The average speed v « 1.22 
is visualized by the thick dotted line. The average flux is exactly at q = 
1 — p — 0.55, as expected from a trivial flow solution in the fundamental 
diagram. Even on the open-road segment the speed does not exceed y™ ax = 
3. Therefore, the bottleneck has no longer an impact on the flow solution. 



on the effective fundamental diagram, which is defined as the average flux 
on the bottleneck loop versus the average density, as shown in Fig.[9l It 
shows three separate values of p (p = 0.0; 0.2; 0.5). To plot the figure, the 
parameters were chosen as L = 1000, Lb = 200, v max 5 and v^ ax 2. 
As the dawdling probability is increased, there is a decline in the flux for 
both parts of the road. This is expected since the vehicles are more likely 
to reduce their speed randomly as p increases. The flow also becomes more 
unstable. 

When p = 0, it can be shown analytically that the effective FD is 
piecewise linear and consists of three segments, joined at the critical values 
Pel — 13/75 w 0.1733 and p C 2 = 1/3. These three linear regimes correspond 
to the three wave patterns found in the previous section. The critical values 
for p are determined in a fashion similar to the critical density values in 
section [3l The linear function 

Q(p) = j^P (ii) 

for < p < 13/75 is also derived from conservation of cars and zero wave 
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speed. The intermediate constant function, determined by the maximum of 
qB for case 1 in Fig. El then joins onto the FD of the open road (no bottleneck; 
Lb — 0). Hence, we have a complete understanding of the effective FD for 
a bottleneck with p = 0. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Density (p) 



Figure 9: Effective fundamental diagrams (dashed lines) for three different 
values of the dawdling probability p from a simulation with L = 1000, 
Lb — 200, v max = 5 and v^ ax — 2. Three distinct regimes can be identified. 
For comparison, the FD for the same simulations without a bottleneck (solid 
line; Lb — 0) exhibit only two distinct regimes. 

However, for p > we resort to numerical results only which again 
exhibit three distinct regimes, or phases, with rather blurry transitions be- 
tween them. The occurrence of such three phases has been discovered in 
previous C A simulations and will we address some of these in section El 

5 Link to optimal- velocity models 

While the above first-order analysis explains the transition between the three 
wave patterns which are observed in the numerical simulations, it does not 
answer why they are selected or why the remaining wave patterns in Fig. [3] 
do not occur. In this section, we will draw a parallel to the microscopic car- 
following optimal- velocity (OV) model. We will simulate a bottleneck in the 
OV model in a way similar to the CA model and analyze how it compares 
to the CA results. It turns out that the same three wave patterns appear 
again and occur exclusively. Using the corresponding continuum model, we 
can explain this phenomenon with second-order analytical methods. 
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5.1 Bottleneck simulation in the OV model 



In the optimal- velocity model [22], the acceleration of each vehicle is de- 
termined by an equation of motion that includes an optimal (or desired) 
velocity V which depends on the distance, h n = Ax n , to the preceding 
vehicle. 

The governing equation reads 

x n = a[V(Ax n ) - x n ] (12) 

where Ax n = x n +i — x n , and n denotes the index of each vehicle at position 
x n . The velocity is represented by x n and acceleration by x n . 

In order to be consistent with the CA model, we use a piecewise linear 
optimal-velocity function with a maximum velocity in the bottleneck at 
60% of the maximum velocity on the open road, again shown qualitatively 
in Fig. [21 For simplicity, we used v™ ax = 0.6 and v max = 1.0 with an offset 
h! — 0.1, representing the distance a vehicle occupies in a jam. Note that 
we are working with a dimensionless model just as it was the case for the 
CA model. 

This results in an optimal-velocity function (bottleneck: subscript B; 
open road: subscript 6) 



0.6 ; h > 0.7 

/i - 0.1 ; h<l.l, 
1.0 ; h > 1.1 



t ,,,n I h — 0.1 ; h < 1.1, , 1 
Vo(h) = { in ' ^ ' (14) 



and the fluxes 



/ l-0.1p ; P > 1/0.7, 
qB = \ 0.6p ; p < 1/0.7 (15) 

a -| 1 -° Ap ; P>1/L1 ' (16) 

q ° \ p ; p< 1/1.1. 1 bj 

The sensitivity a is chosen so that we guarantee linear stability. The flow 
stability is determined by [22] 

V'(/0 < a/2, (17) 

which also defines a range of headways for which the model is unstable. For 
the simulations performed, a = 2.0 was chosen to maintain stable flow. 

Initially, we place N vehicles randomly along a loop of length L with a 
bottleneck of length Lb and let the system evolve. 
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5.2 Wave selection in the fundamental diagram 



Using the above equations for the flux in the bottleneck and on the open 
road, the fundamental diagram is (qualitatively) identical to that of the CA 
model in Fig. [3J After a variety of simulations are carried out, the same 
three wave patterns which occurred in the CA model, also emerge in the OV 
model, while other wave patterns never appear. These three wave patterns 
also form three traffic phases which are qualitatively identical to those of 
the CA model (see Fig.©. 

We find: i) Two shocks at low density, one occurring at either end of 
the bottleneck, corresponding to case 2 in the CA model; ii) two shocks at 
medium density, one occurring at the upstream bottleneck boundary and 
one (classical Lax shock) appearing on the open road (case 1 in the CA 
model); hi) uniform, homogeneous flow (case 3 in the CA model). 

These numerical results, including the critical densities of each regime, 
can again be discussed with first-order analytical methods, namely conser- 
vation of cars and zero wave speed, i.e. kinematic wave theory. This analysis 
explains again the connection between the wave patterns observed in the 
(OV) model but it does not explain why these are the only wave connec- 
tions selected in Fig.[3l In order to investigate this phenomenon in more 
detail, we will now use continuum theory and analyse the travelling-wave 
phase plane of the continuum analogue of the OV model. 

5.3 Travelling- wave phase-plane analysis 

Previous work in continuum theory by Berg et al [12] and Lee et al. [13] 
will now be used to analyse the travelling-wave patterns expected from the 
fundamental diagram of the OV model and, hence, the CA model. A priori, 
it is unclear whether wave patterns may occur in the OV model which cannot 
be found in the CA model. Hence, second-order continuum theory will be 
applied to analyse all plateau connections in the OV model. 

Lee et al. [13] derived a second-order continuum model of the OV model, 
using density p and velocity v as the two variables. Using the conservation 
of cars 



Pt + (vp) x = 



(18) 



and the dynamic equation involving the optimal velocity 



vt + vv x = a[V(p) -v} + 



aV'(p) 
2p 




a 



2 u xxi 



(19) 
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the stability of the shock tails will now be analysed, meaning the stability 
of the equilibrium points (EP) in the corresponding travelling-wave phase 
plane. Here, we define V(p) = V(l/p) = V(h). 



5.3.1 Phase-plane analysis under constant optimal velocity 

A constant optimal velocity, V'{p) — 0, refers to the positive slopes of the 
fundamental diagram where the flux is found from q = y max p in Eqs. ffT5|) 
([TOD . Therefore, we have V(h) = V(l/p) = V(p) = v max and V'(p) = 0, 
where v max is the constant maximum optimal velocity. Note that the same 
analysis is applied to the bottleneck regime by substituting v^ ax for y max . 
With these parameters and assuming stationary shock waves (v t — pt — 
0), Eq. (H9j) becomes 



vv x = a[v™-v] + -^v xx . (20) 

Using the definition for the flux, q — vp^ and linearization via v = v max + y 
yields 

ni .max2 

v ™*v x = -av + -^v xx . (21) 

Here, the flux q is given by the boundary conditions at infinity. 

This equation can now be written as two first-order ordinary differential 
equations 

v x = w, (22) 
^ = J^2^- + H, (23) 

where we have dropped the "hat" notation for v. 

To find the stability of the equilibrium points, the eigenvalues of the 
preceding equations are found by rewriting the equations in matrix form 





6q 2 




(24) 



Trying the ansatz v,w ~ exp(Ax), this yields the following two roots of the 
characteristic equation 



A 12 = J2L± JfJf-V + Jf,. (25) 

CLV max V V CLV max J ymaxZ v / 
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Since g, v max > holds, the root is always larger than the term preceding 
the root and the eigenvalues represent a saddle in the phase plane of v and 
v x . This allows wave connections both i) from an upstream plateau (equilib- 
rium point) to this equilibrium point, representing the second plateau, along 
the stable manifold and ii) from this plateau to another plateau along the 
unstable manifold. Hence, wave connections are possible along the positive 
slopes of the fundamental diagram, as in case 2 in Fig. El which can now be 
interpreted as a saddle-saddle connection. 



5.3.2 Phase-plane analysis for an optimal velocity dependent upon 
headway 

An analysis is now performed for the monotonically decreasing part of the 
fundamental diagram where the OV function does vary with headway and, 
consequently, with density. Here, we have V{h) = V(p) = l/p — 0.1. The 
same analysis is carried out as before, however, this time for the dynamical 
equation 



VVn, 



- -0.1 -v 
P 



2p3 Px + 6pi Vxx ' 



(26) 



Using q — vp again, and hence p x — — p^-, we obtain 



vv x 



av av 
+ 2^ Vx + 6^ Vxx ' 



(27) 



Now let v = v + v and v — y q _i ? we find 



. av ^ av ^ 

1)V + TT^V X + —nV x 



2q 2 



(28) 



Writing this again in terms of two first-order differential equations, we obtain 



w x 



6q 2 



v 1 



2g 2 



w — a [ 1 } v 



(29) 
(30) 



To find the stability of the equilibrium points, the eigenvalues of the pre- 
ceding equations are determined via matrix notation 







V 

w 



(31) 
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and so 




(32) 



Since < q < 1 holds, it can be seen that ^1 — ^ < and the argument 
under the square root can be positive or negative. In any case, stability 
is governed by the sign of the real part of Xi^ and that is always deter- 
mined by ^1 — T^p-^ since the root is either purely imaginary or real but 

smaller than ^1 — ^p-)- Since q,a,v > holds, the real part will either 

be greater than or less than zero depending on the term ^1 — ^2 ) • If a > 2 
and < q < 1, which is the case in our simulations (see Eqs. (fT5l) and (TT6l) 
where q < 1/1.1), the real part will be less than zero. Hence, the equilibrium 
point (EP) is a stable node (or spiral). Therefore, we can connect from an 
upstream plateau to this EP but we cannot connect from this EP to an- 
other equilibrium point downstream to form a wave pattern. This prevents 
the formation of patterns 3, 4 and 5 in Fig. El which would involve a non- 
permissible connection starting from an EP of the monotonically decreasing 
branch of the FD. 

For the same reasons, case 1 in Fig. [3] must be considered as a limiting 
scenario of the above-mentioned saddle-saddle connection of case 2 between 
the monotonically increasing branch of the bottleneck FD and the open-road 
FD, respectively. 

Future work will investigate a bottleneck simulation for unstable flow 
where a < 2 holds. Then it is possible to find unstable nodes in Eq. (l32l) 
since we can have 



i -2? >i 

or, equivalent ly, 



(33) 



a < 2q z 



(34) 



5.3.3 Further thoughts 

At this stage, it is unclear how the instability of the flow will affect the wave 
patterns and whether some stable wave patterns will persist. 
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Also, it should be noted that a similar bottleneck simulation with the 
(linearly stable) OV model [22] has been carried out by Wilson et al. [2^ (25]. 
Here, the bottleneck was simulated by a reduction factor in the optimal- 
velocity function, resulting in new wave patterns. A coarse-graining method 
reveals that second-order theory might be insufficient to explain all wave 
patterns at the microscopic level. However, a higher-order continuum model 
is not available. Moreover, if one follows the idea of Berg et al. [12] and 
extended their continuum model to higher order, it would involve second- 
order (and higher) derivatives of V(p) and they all vanish for a piecewise 
linear OV function so that a higher-order modelling approach would fail in 
the context of this paper. 

Furthermore, the continuum model of Berg et al. completely fails in ex- 
plaining the saddle-saddle connections since one ends up with an insufficient 
first-order model 

w x = a[v max - v] (35) 

along those branches of the FD. This clearly shows the advantage of the Lee 
et al model, adding to its advantage over the Berg et al model in terms 
of matching the linear stability criterion of the discrete OV model at short 
wave lengths. 

In summary, kinematic wave theory alone cannot fully explain the se- 
lection of stationary shocks in the bottleneck simulation. Surely, stationary 
wave patterns can be expected at the boundaries due to a non-smooth change 
in model parameters. However, it is the dynamics of the model which deter- 
mine the actual observable patterns, and these resemble each other in the 
OV and CA model simulations. This includes the appearance of three traffic 
phases in either model. 

6 Comparison to open-system and other 
non-homogeneous CA models 

There have been several publications about CA models in the literature 
which contain either two [7J , three [HI El El El [9] or several traffic phases [4J , 
as defined in section [H 

Yukawa et al [3] simulated a blockage on a loop in the rule-184 CA 
model, having a maximum speed y max = 1. At exactly one site along 
the loop, the hopping probability of an otherwise deterministic model is 
reduced from r = 1 to < r < 1. The authors find three traffic phases, 
depending on the average traffic density and the hopping parameter r. In 
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order of increasing density, they are: i) free, ii) constant-flow and iii) jam 
phase. This resembles the three wave patterns in this work (see section [4]), 
where the constant-flow phase is determined by q^ ax = 3/4 in the regime 
0.17 < p < 0.25. However, while the free phase in the rule-184 CA model 
is the same with or without a blockage site, this is not the case in our 
simulations. In addition, the present work considers a fully deterministic CA 
model and focuses on its emerging wave patterns. Moreover, the bottleneck 
simulations in section [3] will change as Lb — ► 1 due to finite size effects since 
the formation of plateaus within the bottleneck will disappear. Therefore, 
planned future work will investigate how the limit Lb — ► 1 will affect the 
results both for y™ ax = 3 and y™ ax = v max — 5 5 an( j a non-zero dawdling 
probability p > in the bottleneck only. 

Lakatos et al [9J extended the idea of Yukawa et al to multiple blockage 
sites in a totally asymmetric simple exclusion process (TASEP) model on 
an open road. Approximate mean-field theory provides good agreement be- 
tween the numerical results and Monte Carlo simulations. Phase diagrams 
are derived, which contain again three phases, but wave patterns are not 
discussed due to the lack thereof. More importantly, the boundary condi- 
tions mainly drive the dynamics of this open system, which brings us to the 
next set of simulations. 

The stationary wave patterns of the present bottleneck simulation would 
change dramatically if the loop was replaced by an open road with fixed but 
different up- and downstream boundary conditions [H El El El E]- The con- 
straint of conservation of cars, which enables the analytical computations of 
phase boundaries between the three wave patterns (see also [3]), disappears 
and the analysis of the model becomes more challenging. Depending on the 
model, three [5j or more phases [4J can be found. Here, stationary wave 
patterns are determined by the boundary conditions alone and can even be 
triggered in an otherwise stable model, with or without a bottleneck. In 
contrast, such patterns may or may not exist on the loop, depending on the 
stability of the model. This fundamental difference is exhibited for the OV 
model on a loop and open road by the work of Bando et al [22] and Berg 
et al. [26], respectively. Therefore, future work will also investigate the bot- 
tleneck setup of the present paper under open-system boundary conditions 
similar to 01 El El El El . 
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7 Conclusion and future work 

In this paper, the dynamics of a bottleneck simulation exhibit a link between 
cellular automata and optimal-velocity traffic models. Three wave patterns, 
which correspond to three distinct traffic phases in the fundamental diagram, 
occur in either model and are qualitatively identical. This is surprising in 
some sense since it connects a model, which is discrete in space and time, to 
a model, which is continuous in space and time. The emerging wave patterns 
were analysed and explained with second-order continuum theory. 

Future work will focus on a CA model with non-zero dawdling probabil- 
ities < p < 1. However, this entails the formation of jams, which overlap 
with the stationary wave patterns. On the other side, similarities to the 
work of Yukawa et al [3] can be expected to emerge from these simulations. 

Moreover, the stationary wave patterns which occur on a long straight 
road, containing a bottleneck and different up- and downstream boundary 
conditions, will be studied. Here, travelling waves cannot interact multiple 
times with stationary patterns. This should give rise to new wave solutions 
as compared to the loop, determined by the boundary conditions. This 
resembles research presented in [5j O [8] but is fundamentally different from 
what was presented in this work. 
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